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Based on state-of-the-art density-functional-theory methods we calculate the stacking-fault energy 
of the paramagnetic random Fe-22.5at.%Mn alloy between 300-800 K. We estimate magnetic thermal 
excitations by considering longitudinal spin-fluctuations. Our results demonstrate that the interplay 
between the magnetic excitations and the thermal lattice expansion is the main factor determining 
the anti-Invar effect, the hep-fee transformation temperature, and the stacking-fault energy, which 
is in excellent agreement with measurements. 
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The intrinsic stacking-fault (SF) is one of the simplest 
planar defects of the crystal lattice. It is characterized 
by a fault in the usual ABC stacking sequence of the 
fee structure, ...ABCAB\ABC..., which resembles locally 
the stacking sequence of the hep structure. The energy to 
create a stacking-fault, the stacking-fault energy (SFE), 
is related to the ductility of the material, but the connec- 
tion between the SFE and the plastic deformation mech- 
anism has so far only been recognized empirically, for 
instance, in high-Mn steels [H, Nevertheless, the exis- 
tence of such a connection provides a unique opportunity 
for a first-principles theory to be a part of intelligent de- 
sign of new materials, and in particular high-performance 
steels, since the SFE is an atomic-scale property readily 
accessible via ab initio methods. 

The ab initio calculation of the SFE in real alloys at fi- 
nite temperature presents two main complications: First, 
the proper description of the atomic-scale structure of 
multicomponcnt alloys is complicated because the atomic 
configuration of alloy components is usually unknown, 
and, moreover, it is still quite cumbersome to model. Sec- 
ond, the accurate treatment of the thermal excitations 
of lattice vibrations and magnetic fluctuations in alloys 
still remains a challenge for any theoretical approach. 
The magnetic fluctuations are particularly important in 
fee and hep Fe-based alloys, which are weakly itinerant 
magnets. In particular, this concerns the technologically 
important Fe-(20-30) at.%Mn alloys, which exhibit both 
Invar and anti-Invar effects, as well as martensitic trans- 
formations and shape memory effects 0, Q ■ 

Special attention is given to the fee Fe-22.5 at.% Mn, 
which exhibits the lowest room-temperature SFE of the 
fee Fe-Mn alloys [H, and thus is of particular interest 
for the automotive industry. This alloy undergoes a 



structural transformation between 7-austenite (fee) and 
e-martensite (hep) which, on cooling, is characterized by 
the martensite start temperature, M s = 375 K Q and, 
on heating, by the austenite start temperature, A s = 
450 K Q. Important is that the Neel temperatures of 
both fee (360 K 0) and hep (230 K 0]) are below M s 
and A s , which implies that the fec-hep transformation 
occurs in the paramagnetic state. In addition, the anti- 
Invar effect, or enhanced thermal lattice expansion, that 
has been correlated to the magneto- volume coupling 0], 
also takes place in the paramagnetic state. It is, thus, 
clear that methods capable of treating the paramagnetic 
state in a sufficiently reliable fashion are required for a 
proper description of any finite-temperature properties 
in this system. 

In this Letter, we study the SFE of the 
Fe-22.5 at.% Mn binary alloy between 300-800 K, 
using a simple but powerful model based on density- 
functional theory to take into account spin fluctuations 
of Fe and Mn atoms, which allows us to obtain excellent 
agreement with experimental data [f|. The key feature 
of the methodology used is the capability to describe the 
strong magneto-volume coupling at finite temperature. 
We show then that the magneto-volume coupling is 
responsible for the anti-Invar effect and the temperature 
dependence of the SFE in the Fe-22.5at.%Mn alloy. 

Two recent ab initio calculations of the SFE in FeMn 
alloys [H,[i|, that were done at zero-K, failed to reproduce 
the finite-temperature experimental values even qualita- 
tively. The main motivation of this work is to show that 
the accurate treatment of the thermal excitations of the 
magneto-volume coupling is essential for the calculation 
of the temperature-dependent SFE for the fee Fe-based 
alloys, in particular Fe-22.5at.%Mn. 
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The SFE is evaluated using the axial next-nearest 
neighbor Ising (ANNNI) model [nj: SFE(T) = 
G hcp (T) + 2G dhcp (T) - 3G fcc {T), where G(T) is the 
Helmholtz free energy, and the volume per atom of the 
ideal hep and dhep structures is the same as for the fee 
structure, implying that only the fee volume needs to be 
known. The alloy configuration is assumed to be com- 
pletely random, which is supported by recent experimen- 
tal data [m j an d the Mn concentration at the SF is the 
same as in the bulk, i.e., the Suzuki effect fl2| is inoper- 
ative due to the very slow diffusion of manganese in the 
austenite phase [l3j . The ANNNI model is known to be 
less accurate than direct supercell calculations, where the 
SF structure is treated explicitly and the geometry can 
be relaxed [HJ]. However, accurate SFE supercell calcu- 
lations for random alloys in the paramagnetic state are 
currently impractical. 

Generally, the free energy is a difficult quantity to eval- 
uate. Using the coarse graining of the partition func- 
tion, the contributions from electronic, magnetic, vibra- 
tional and configurational excitations can be separated 
out based on their different time scales. Clearly, the 
entropy contribution from magnetic excitations in the 
paramagnetic state can be significant and needs to be as- 
sessed. As for the vibrational entropy, there are currently 
no available theoretical tools to determine it accurately 
in paramagnetic random alloys. However, its contribu- 
tion to the free energy differences required to evaluate 
the SFE has been argued to be small [Hj]. The config- 
urational entropy does not contribute to the SFE in the 
case of a completely random alloy. Lastly, we regard the 
effect of the local lattice relaxations to be negligible since 
the atomic radii of iron and manganese are similar. 

Thus, the magnetic excitations seem to be the main 
entropy contributor to the SFE. The relevant part 
of the Helmholtz free energy per site is then deter- 
mined as G(T) = F a (T) - TS mag (T), where F el {T) = 
F e i \Sws(T) 7 {fhi(T)}} is the electronic part of the free 
energy, which includes one-electron excitations and de- 
pends on the volume given in terms of the Wigner-Seitz 
(WS) radius Sws(T) and the average local magnetic mo- 
ments rhi(T) of each alloy component i; the magnetic 
entropy S mag (T) = S mag [{fhi(T)}] is determined by the 
thermal fluctuations of the magnetic moments. 

We start with a detailed discussion of the equilibrium 
volume at finite temperature Sws (T) and its relation to 
the magnetic moments fhi(T). First, we consider the 
equation of state at T = K of the fee and hep ran- 
dom alloys in the paramagnetic configuration. In order 
to provide an accurate account of the local environment 
effects and energetics, we model the alloy by 384-atom 
supercells consisting of 298 Fe and 86 Mn atoms dis- 
tributed randomly on the underlying fee and hep lat- 
tices, i.e., the alloy composition is Fe .776Mn .224- The 
DFT self-consistent total energy calculations are done in 
the disordered local moment (DLM) [l6| configuration by 
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FIG. 1. (Color online) Total energies (upper-left panel) and 
zero-K DLM magnetic moments (lower-left panel) of Fe and 
Mn in fee (circles) and hep (triangles) Feo.77eMno.224, as func- 
tion of volume. Two dashed vertical lines show the experimen- 
tal WS radii for this alloy at 300 and 800 K Helmholtz 
free energies (upper-right panel) and finite-temperature DLM 
moments (lower-right panel) of the fee Feo.776Mno.224 alloy at 
400 K (squares) and 800 K (diamonds). The total energy of 
the Feo.776Mno.224 alloy at zero-K is shown by filled circles. 



using the locally self-consistent Green's function (LSGF) 



method 121 ll3 implemented within the exact muffin- 



tin orbital method combined with the full-charge density 
technique (EMTO-FCD) 0, called LSGF-EMTO [H. 

In the upper-left and lower-left panels of Fig. [TJ we 
show, respectively, the total energies and the DLM local 
magnetic moments of Fe and Mn as a function of the WS 
radius for the fee and hep Feo.776Mno.224 alloys. We ob- 
tain a zero-K fee equilibrium WS radius, 2.563 a.u., which 
is much smaller than the room-temperature experimen- 
tal value, 2.658 a.u. [23|, indicated by the dashed vertical 
line in Fig. [T] 24|. Our zero-K result is in agreement with 
a recent full-potential supercell calculation, 2.574 a.u. [§], 
while at variance with a previous EMTO-FCD study, 
2.605 a.u. [25)]. The latter calculation, however, used the 
frozen-core approximation, a smaller set of k-points, and 
had a freedom in the choice of unspecified screening pa- 
rameters, which can explain the disagreement. 

From Fig. [2 we realize that both fee and hep phases 
possess a low-spin equilibrium volume and exhibit a low- 
spin (LS)- high-spin (HS) transition with increasing WS 



radii, indicative of the anti- Invar behavior |26[. In ad- 
dition, the fee structure becomes more stable than hep 
above 2.660 a.u., while simultaneously, a LS hcp-HS fee 
transition takes place, demonstrating the coupling be- 
tween the structural transformation and anti-Invar be- 
havior in the paramagnetic state, which has been dis- 



cussed in the literature |26| . 

At finite temperature, the Hclmholtz free energy must 
be calculated instead of the total energy. The electronic 
single-particle excitations are described via the Fermi 
function, while we use the simplest model to describe the 
longitudinal spin-fluctuations (LSP) in the DLM state. 
Specifically, we prescribe to each site p in the super- 
cell, having DLM magnetic moment m p , the entropy 
Smag = fcs hi(m p + 1). The value of m p is determined 
by the minimization of the Hclmholtz free energy, and in 
this way, it corresponds to the average magnitude of the 
magnetic excitations in the DLM state. We note that a 
more elaborate approach to the treatment of LSF will be 
used below for the calculation of the SFE. 

We have calculated the free energy as a function of the 
WS radius at 400 and 800 K using LSGF-EMTO, taking 
into account the LSF as described above. Compared with 
the zero-K result, the minima of the free energy curves 
arc significantly shifted towards larger volumes upper- 
right panel in Fig. [1]), which is caused by the drastic 
increase of the magnetic moments due to the LSF (lower- 
right panel in Fig. [1]) . From these data, one can estimate 
Sws(T) by taking into account the vibrational entropy 
using the Debye-Gruneisen (DG) model |27|. We obtain 
2.650 and 2.698 a.u. at 400 and 800 K, respectively, in 
agreement with the corresponding experimental Sws(T), 
2.662 and 2.682 a.u., estimated from X-ray diffraction 
(XRD) [23| and dilatometric measurements @. 

We also obtain, however, a very large isotropic 
Griineisen constant, about 4, and a strongly overesti- 
mated value of the thermal-expansion coefficient, which 
shows that the DG model cannot be used for the quanti- 
tative determination of Sws{T) in systems with a large 
magnetically-induced anharmoncity. Nevertheless, the 
model shows the strong magneto- volume coupling in this 
system and allows us to capture qualitatively the crucial 
feature of the anti- Invar effect, namely, enhanced ther- 
mal lattice expansion caused by the so-called moment- 
volume instabilities [26|, i.e., the thermally activated LS 
low-volumc-HS high-volume transition. 

The results show that the magneto- volume coupling is 
an essential feature of the FeMn alloy, and it is, thus, 
important to have accurate values of Sws{T), and the 
average local magnetic moment of each alloy component 
i, fhi(T) = rhi[T; Sws(T)], at a given temperature T, for 
a proper description of the system. There are no theo- 
retical methods available to accurately calculate Sws(T) 
for paramagnetic random alloys, but it can be measured. 
The situation is, however, the opposite with respect to 
the magnetic moments fhi (T) in the paramagnetic state, 
which are difficult to measure, 

but can be calculated by the procedure described be- 
low. We, therefore, combine the values of Sws(T) es- 
timated from the XRD [i^l and dilatometric measure- 
ments Q with the theoretically evaluated values of fhi (T) 
to obtain the free energy and finally the SFE at finite 
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FIG. 2. (Color online) The energy of the LSF, J;, ob- 
tained for different structures, at the fee Sws(T = 360 K) 
= 2.660 a.u. in the EMTO-FCD calculations. 



temperatures. 

To this end, we proceed to a more elaborate treatment 
of the LSF and determine fhi(T) following the approach 
in Ref. [1^ generalized to the case of a binary alloy. In 
this procedure, fhi(T) are evaluated as the thermody- 
namic averages of the fluctuating magnetic moments of 
iron and manganese at a corresponding temperature. For 
that purpose we assume that the free energy is dominated 
by single-site magnetic fluctuations in the paramagnetic 
DLM state. In this Hamiltonian describing the 

LSF can be defined as H rnag = ^ Ji(rrii)ci, with i =Fe, 
Mn, where Ji(uii) is the energy necessary to excite the 
magnetic moment mi of the ith alloy component. In cal- 
culating JFe(Mn) ; we nx one of the magnetic moments, 
nT'Fe(Mn); while letting another, vn,Mn{Fe)i relax. 

In Fig. [2] we show the results for JArrii) obtained for 
fee S ws iT = 360 K) = 2.660 a.u. 0, in the EMTO- 
FCD calculations |29j. In all the cases, Ji have very shal- 
low minima indicating low-energetic accessibility of spin 
fluctuations, with large variety of the magnetic moments, 
which reflects the weak itinerant nature of magnetism in 
this system. The minima of the Ji(rrii) define the equi- 
librium magnetic moments in the DLM ground state, 



,DLM 



rriiiT = 0), and we obtain (mp 



DLM 



DLM 
l Mn 



= (1.4,1. for the fee, (0.0,0.0) fis for the hep, and 
(1.2,0.8) hb for the dhep structures, respectively. Note 



that these ■ 



,DLM 



are close to those in the lower-left panel 



of Fig. [TJ at Sws = 2.660 a.u., that have been obtained 
at zero-K in the LSGF-EMTO calculations. 

At finite temperature T, the average magnetic mo- 
ments, fhi(T), are readily obtained from the Ji by a 
simple Monte Carlo integration for each component. In 
doing so, we assume the complete coupling between Ion- 
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FIG. 3. (Color online) Average magnetic moments of Fe and 
Mn in different structures (upper panel), and the SFE and 
its components (lower panel) as a function of temperature. 
Vertical line indicates the theoretical estimate of the hcp-fcc 
phase transition. 



gitudinaJ and transverse fluctuations, i.e., the magnitude 
of the local magnetic moment is determined by its x, 
y, and z-projections, when integrating over the classical 
spin phase space. The results are shown in the upper 
panel of Fig. [3] for different structures in the tempera- 
ture interval of 300-800 K, and were obtained using the 
Ji of Fig. [21 Note that our room-temperature rhi(T) 
for the fee phase is in agreement with neutron diffraction 
room-temperature measurements of the average moment 



per site in fee Fe-30at.%Mn, 1.5 [is [31 



Finally, in order to obtain the SFE, we use the calcu- 
lated rhi(T) of Fig. [3] together with the experimental fee 
Sws(T) in the free energy calculations of the fee, hep and 
dhep structures. In the lower panel of Fig.[3]one can read- 
ily see that despite a small quantitative discrepancy with 
the experimental data the temperature dependence 
of the SFE is nicely reproduced. Moreover, the entropic, 
TAS m ag (T) , and electronic free energy, AF e i (T) , contri- 
butions behave quite differently with temperature, with 
the latter being fully responsible for the temperature de- 
pendence of the SFE. The entropic contribution, on the 
other hand, just shifts (although quite substantially) the 



value of the SFE towards the experimental value. This 
result is contrasted with earlier results by Vitos et al. for 
FeCrNi alloys [l5[ , where the temperature dependence of 
the SFE was found to be determined solely by the mag- 
netic excitations. 

A detailed scrutiny of the sensitivity of the SFE to 
the values of Sws(T) and rrn{T) enables us to establish 
the crucial role that thermal excitations and magneto- 
volume coupling play in the behavior of the SFE. In par- 
ticular, neglecting the temperature dependence of either 
the volume, Sws, or the magnetic moments results, re- 
spectively, in a largely underestimated value of the SFE, 
-250 mJ/m 2 at 300 K, when using the equilibrium the- 
oretical zero-K Sws and rrii(T) from Fig. [3[ or in an 
overestimated value of 150 mJ/m 2 at 300 K with the ex- 
perimental Sws(T) and zero-K mf LM . 

The obtained free energies allow us also to investigate 
the structural stability by examining the so-called driv- 
ing force AG(T) = G hcp {T) - G fcc {T). Our value of 
the transition temperature, determined by the condition 
AG (To) = 0, is 375 K, which is obtained using the ideal 
c/o ratio and the same atomic volume for the hep phase 
as for the fee one. It is in reasonable agreement with the 
experimental estimate T cxp w (M s + A s )/2 = 412 K [f|. 

We would like also to remark on the importance of 
including the dhep term to the expression for the SFE. In 
the axial nearest-neighbor Ising (ANNI) model, the SFE 
is equal to 2AG. The comparison of the SFE calculated 
using the next- nearest- neighbor model and 2 AG in Fig. [3] 
makes clear that the contribution of the dhep structure 
is important for the quantitative description of the SFE. 

In summary, our results demonstrate that the interplay 
between the LSF magnetic excitations and the thermal 
lattice expansion is the main factor determining the anti- 
Invar effect, the hcp-fcc transformation temperature, and 
the SFE in the paramagnetic random Fc-22.5at.%Mn 
alloy. In principle, this strong sensitivity of the SFE 
with respect to the quite elusive thermal effects of the 
magneto-volume coupling makes theoretical studies in 
magnetic alloys rather difficult. At the same time, we 
have clearly demonstrated that the accurate account of 
important thermal contributions within the DFT ab ini- 
tio approach allows not only to obtain a good quanti- 
tative agreement with experiment but also to identify 
mechanisms dominating the formation and stability of 
the SF at finite temperatures. The formalism presented 
here can be further generalized to the case of austcnitic 
stainless steels, which are many-component systems with 
the presence of strong atomic short-range order and local 
lattice relaxation effects, thereby opening a perspective of 
the intelligent design of new high-performance materials. 
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